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Negative dimensional integration method (NDIM) is revealing itself as a very useful technique 
for computing Feynman integrals, massless and/or massive, covariant and non-covariant alike. Up 
to now, however, the illustrative calculations done using such method are mostly covariant scalar 
integrals, without numerator factors. Here we show how those integrals with tensorial structures can 
also be handled with easiness and in a straightforward manner. However, contrary to the absence of 
significant features in the usual approach, here the NDIM also allows us to come across surprising 
unsuspected bonuses. In this line, we present two alternative ways of working out the integrals 
and illustrate them by taking the easiest Feynman integrals in this category that emerges in the 
computation of a standard one-loop self-energy diagram. One of the novel and as yet unsuspected 
bonus is that there are degeneracies in the way one can express the final result for the referred 
Feynman integral. 



In an effort to make sense of diverging integrals that popped out in the field theoretical approach to transition 
amplitudes, scattering matrices and so forth, physicists introduced and developed the concept of extended dimensions 
0. This can be interpreted in a pragmatic way as a mere artifact to go around a difficult problem. Nonetheless 
though it may be so, the principle of analytic continuation behind it is mathematically well-founded and established, 
being self-consistent and well-defined. Thus, if we are allowed to say that the beauty and the power of Mathematics 
reside in the possibility of defining abstract entities that have no real connection to our physical world, and from 
such entities which one could rightfully call "outside of this world" we can draw either some sense from it or even 
pertinent and meaningful properties that become relevant to our dimensionality, then the effort to go to these frontiers 
is worthwhile and enriching. In a very personal way of viewing things, we think of NDIM in such terms. 

We then play with negative dimensions and with precise analytic continuations, so that interesting results do emerge 
from our exploration of this kind of realm. A very useful technique stemming from such incursion is the method of 
integrating Feynman integrals in negative dimensions Q . Instead of the usual field propagators in the denominator 
of the integrands, here we have them as numerators. In other words, in essence what we have here are integrands of 
polynomial type. Of course, once the integral is performed in negative dimensions, it must be analytically continued 
back to our real, positive dimensional world. The basis for doing this is set forth in our previous papers |j|-|5| ■ 

Our aim in this work is to further illustrate the methodology of NDIM and for this purpose we take examples from 
the one- loop vacuum polarization tensor diagram which generates some Feynman integrals with tensorial structures. 
We show that the calculation of integrals with tensorial structures can be dealt with propriety using the NDIM 
technology. Moreover, we show that this can be approached in at least two ways, which we consider with details in 
the next sections. A first approach is to just "copy" the steps used in the traditional positive dimensional approach, i.e., 
using derivative identities in the integrands. A second, novel approach is to define the relevant negative dimensional 
integral corresponding to the Feynman integral we want to evaluate right from the beginning and proceed from there. 
Just to make the illustrations simpler and clearer, we restrict ourselves to massless fields, but generalization to massive 
ones is not difficult to do. 
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I. INTRODUCTION 



II. USING DIFFERENTIAL IDENTITIES 



Let us first consider the following (vectorial) Feynman integral: 




(1) 
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which clearly emerges in the calculation of vacuum polarization tensor of, e.g., quantum electrodynamics. This, of 
course, is easily calculated in the standard procedure of positive dimensions. How is this done in the NDIM context? 

The structure of the above integral immediately suggests that a possible way of starting off NDIM calculation is to 
consider a Gaussian-like integral of the type 



Q» = I d 2D q <f e- aq -/>(«-*) . (2) 
This, in terms of the negative dimensional integral j, is therefore given by 

= jt(-iy +1 ^-^(hJ,D;p). (3) 

ij=0 ' 3 ' 

On the other hand, performing the momentum integration of equation (Q) through the use of the following identity, 

^e* P 4^ e * P - ( 4 ) 
2/3 dp^ 



we get 



= p»n D (-ir(-x-D-l)\ ULL Sa+b , x+D+1 , (5) 

x. a. 6—0 



a x+a px+b+1 (p%y 

= yn ^ \~ L ) \—* — u — i.y. 

x,a,b=0 

where A = a + 0. 

Comparison between equations (Q) and (|5|) term by term yields the result for j, D;p). After analytic continu- 
ation to positive dimensions and negative values of exponents (i, j) 0, we get 

= ^^nrr^- (6) 

\—a\2a + D + 1) 

where we have used a = i + j + D and the definition for Pochhammer's symbols 

Hk) - H), - T -^ 1 ■ (7) 

T{a) 



Next, we consider the tensorial Feynman integral 



i2u) 1 1 



q 2 (q - p) 

The procedure is completely analogous, now starting from 



^ = / d ^ -rrrh^ . (8) 



= J d 2D q q» q v e-^-Pii-P) 2 . (9) 



Here we are going to quote only the final result, which reads 



W> \ p p (- a \2a + D + 2) 2 (-a- l\2a + D + 3) ? 
In a similar manner, we can evaluate the following integrals very easily: 
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jfiup _ -jP^P 



[P ' \ 2 (-a- l\2a + D + A) ^ P 1 P (-a\2a + D + 3) J 
where T pvp = p p g vp + p v + pP g>* v , and 



inept, _ -rHupt; 
1 ~ ^AC 



where 



4 Afivpi 2 K>pup<; 

— A — + l —— — v B + pp p ypp" t p 



r B = 
v v = 



(-i\a + 2) (- 




-2) 


(-<T- 2\2a- 




6) 


(-i\a + l)(- 


-ik-t 


-3) 


(-<7- l\2a- 




5) 


Hk) (-ik 


+ 4) 




(-ct|2o- + D 


+ 4) 





(11) 



(12) 



(13) 



with A^ VIK = g plJ g pq + g w g vq + g>* g vp and B pv(K > = p p p v g p<i + permutations. All these results agree with those 
in the Apendix A of M. 



III. USING PURE NDIM TECHNIQUE 

In order to calculate tensorial structures in Feynman integrals, we can adopt another alternative approach. Let us 
consider the following integral: 

J= [d 2D q P/^' ;>0 (14) 



V {q-p) 2 

Of course, for I > the tensorial structure is implicit, being contracted with external vector p. The advantage of 
this approach is that it takes care of all the tensorial structures at the same time. 

So, instead of using equation (Q) or equation ([)]), etc. as our starting point, the structure of the Feynman integrals 
in equation ([[]) and equation (||) also suggests another possible way of defining the Gaussian-like integral of interest 
to begin with, namely, 

H = [ d^qe-^-^-rf-^ 2 ^ . (15) 



This defines the negative dimensional integral I, D;p) as follows: 

, y Ri V 

i,j,l=0 



n= f; fd 2D q ( q 2 y[( q - P ) 2 Y(2 q - P ) 1 



= E ^r +3+l ^fj^J^D;p). (16) 

i,j,l=0 ' J ' 

On the other hand, from ( |l5| ) we get also 

H = n D jr (-l) x+y 2y ( ~ g '^f^ a x+a (3 x+ y+ b 1 y+ 2z (17) 

x,...,b=0 
a+b=-a' -D 

where a' = x + y + z = i + j + I + D, or a' = a + /. 

Therefore, the solution for J(i, j,/, D;p) is obtained from the solving of a system of linear algebraic equations of 
the following form j| : 
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I = x + a 

j = x+y+b 

I = y + 2z 

a' = x + y + z 



(18) 



It is very easy to see that the above system is formed by four equations but five "unknowns" (the sum indices 
x, y, z, a, b). Therefore, it can only be solved in terms of one of the sum indices x, y, z, a, or b. For each of these 
remnant indices, the sum yields 3F2 hypergeometric functions of unit argument, as follows: 



J { A S F } =A {S} 3 F 2 (a, b, c;c, f|l) 



where the set {S} = {x, y cvcn , y a dd, z, a, with 



{-j\2i + 2j + 1 + 2D) 



(i + j + D\ I + D) (1 + l\ 2i + 2j + 2D) ' 



(19) 



(20) 



(-i\2i+y + D)(-j\2j+il + D) 



A ovcn = ir D (v 2 Y 

v V ' {-i-j-\l-D\2i + 23 + ll + ZD){l + l\-\iy 



(21) 



A odd = _ 2 A cven + d \ 5) H -J-¥-D\ |) |) 



(22) 



A*=tt d (2p 2 )V) i+i+B 



(-i|2i + / + £>)(-j|2j + £>) 
H-.7-D|2i + 2j + Z + 3L>)' 



(23) 



A a =7r c (2p 2 ) i (p 2 )^+ I) (-4)^ 



c (-j|2j + J) 
(1 + I|2j + 2D)' 



(24) 



k d {p 2 Y' (-1)' (-i| -i-j -l-2D)(-j\2i + j + 1 + 2D) 



A h = 



2^+ l + 2D (l + l\t + D) 

and the corresponding parameters of hypergeometric functions given by: 



Parameters 


3^1 


pa, oven 


3^2 


a 


—i 


i + ±l + D 


i + ±l + D + ± 


b 


-i-j-±l-D 


-i-j-M-D 


-i-j-±l-D + ± 


c 


-i-j-tf-D + i 


—if 

2 l 


--1 + - 


e 


1-i-j-D 


1-j-U-D 

—T 2 




f 


1 - 2% - j - I - 2D 


2 


2 



Parameters 




3^1 




a 


j + D 


—i 


i + j + I + 2D 


b 


-¥ 


i+j+l+2D 


i + ±l + D 


c 


-v+i 


J+D 


i+±l + D + ± 


e 


l + i + j + D 


: j + ±l + D+± 


1 + i + l + D 


f 


1-i-l-D 


l+j + y + D 


1 + 2i + j + I + 2D 



(25) 



1 Note that the y index has been split into its even and odd sectors. 
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Observe that in the process of analytic continuation to our physical world (D > 0), exponents i, j are analytically 
continued to allow for negative values, whereas the exponent I must be left untouched, since, by definition I > in 
the original Feynman integral . 

One of the interesting features of the NDIM technique is that it can give rise to degenerate solutions for the same 
Feynman integral. All the answers we have above, although seemingly distinct, are in fact only different ways of 
expressing the same thing. This means that for this particular case, where the solutions are degenerate, taking one of 
them will suffice. The equivalence of the different forms in which the solutions are expressed is shown in the appendix. 

Since we have this freedom of choice, looking at the hypergeometric functions whose parameters are listed in the 
table, we see that the most convenient solution is given by the one coming from solving the system in terms of the 
summation index z. The reason for this is based on the fact that two of its numerator parameters, namely b = — hi 
and c = — hi + i, readily leads to truncated series for I = even and / = odd respectively. Then 

J = 

= n D (pY 2 l H|2i + 1 + D) (_i|2j + D) 



{-i- j - D\2i + 2j + 1 + 3D) 
x 3 F 2 (j + D, ~l, + l-i-l-D, l + i + j + D l^j . (26) 

Now it remains for us to check the results we obtained so far by assigning explicit values for the exponents i, j, and 
I in (^), (|l0|), (|Tl|), ( |l2| ) and (|2^). Let us begin with i = j = — 1 in @, and in order to facilitate the comparison, we 
shall compute 

2„n-i,-i, i >, f )^.V)°- r(p)r ff 2 -!^ (2 - D) ■ <m 

This result is to be compared with the one coming from equation ( p6[ ) for the particular case when i = j = — 1, 
and 1 = 1. It can be seen straight away that the numerator parameter — -A + , of the hypergeometric function 3F2 
vanishes for the particular value I = 1, so that only the first term, that is one, in the series defining it is relevant, and 

J(-l, -1, 1, D; P ) = 2 (pT- 1 (1 M (1| ^ 2) , (28) 
K , , , ,yj \f ) ( - 2 _ D | 3D _ 3 ) ' v I 



(29) 



which is, of course, exactly equal to equation ( |27| ) as it should be. 
From (O) we get (after contracting it with Ap^p^): 



v ' ' ,y> KF ' \ 2D) T(2D) 



This now is to be compared to the result coming from equation ( |26|) for the particular case when i = j = — 1, and 
1 = 2. 



J(-l, -1,2, D;p) =4n° ( P T ^ OJg Z j) ^ U -- , -D 
y ' ' ' ,yj yF ' (2 - D\3D -2) \ ' 2 



(2-D\3D-2) { 2D 

Note that the a and /parameters coalesce into the same value D — 1, so that the 3-F2 becomes a 2-F1 hypergeometric 
function. Moreover the numerator parameter b = —hi turns out to be a negative integer unity for I = 2, so that the 
series truncate at the second term, and the final result is exactly equal to the RHS of equation (p9[). 

In a completely analogous way, we get from ([ll]) contracted with 8 p^PvPp 

while from (|2^) with i = j = — 1 and I = 3 we get 

j ( -i,-i,3,Ap)^8^( P T +l(1 | J + 1 ? (1|J " 2) ^f--, -1; -D-i 

yF ' T(2D + 1) \ 2(D + 1) < 
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which is exactly the same as ( |3~l| ) above. 

Finally, from ([l2] ) contracted with Wp^p^ppp^ we get 



x 1- 



D + 2 4(D + !)(£> + 2) 



while from rt2q) with i = j = — 1 and I = 4 we get 



,(-!, -1,4, - eiywg,- V , (-,, 4 - 

= 1fi7r ^^ +2 (l|-P + 2)(l|-P-2) f _ _J_ 3 

^ J (2-£>|3£>) \ £> + 2 4(£) + !)(£> + 2) 



(33) 



(34) 



in complete agreement with (p3) above. 



IV. CONCLUSION 



We have shown in this paper how we can work out Feynman integrals with tensorial structures in the context of 
NDIM. There are two equivalent approaches for doing this, either by using differential identities to work them out 
one by one (vector, rank two tensor, and so on) mirroring the positive dimensional technique, or by using pure NDIM 
methodology to get simultaneous results. The former technique does not bring any new feature while the latter one 
allows us to get this new feature of degenerate solutions, plus the bonus of having them all at once. As we have 
noticed before [Q, the pure NDIM methodology shows itself more powerful in that it gives equivalent forms of a 
six-fold degenerate solution for the integral. And not only this, the solutions we get are simultaneously obtained. 
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V. APPENDIX 



In this Appendix, we shall show in great detail the equivalence of the six solutions generated by the pure NDIM 
technology in the computation of the one-loop Feynman integrals with tensorial structures in the numerator. In order 
to do this, let us first consider the solution Jf c ■ Its corresponding ^F^ hypergeometric function has parameters 
given in the table of section III. It is clear from its parameter b z = —^l that for I = even = 2m, m = 0, 1,2, ... 
the hypergeometric function is actually a truncated series. In a similar way, from its parameter c z = — -A + h for 
/ = odd = 2m +1, m = 0, 1, 2, ... the hypergeometric function is also a truncated series. 

Truncated hypergeometric series of the form P F q , p = q+ 1 can be inverted |^,^| from variable, say, x m to x ~ 1 ■ F° r 
the particular case of 3F2 with unit argument, such inversion leads to an identity between them. This is expressed in 
a shorthand notation as 

r(e - c)r(/ - c)F p (0; 4, 5) = (-)"T(1 - a)T(l - 6)F„(3; 1, 2) (35) 

where 

F P (0;4,5)= 3fa ^ C .jw^ |1) . (36) 



and 



r( s )r(e)r(/) ' 



3-^2(1 + c- e. 1 + c-f, c: 1-a + c, 1 — 6 + c 1 1) , , 

1 (s)I (1 — a + cjl (1 — + c) 



G 



with to denoting the negative integer numerator parameter, and s = e + f — a — b — c. 
Another way of writing this up is (for general variable x) 



3 F 2 (-m, a x , a 2 ;c, px\x) = ®m^F 2 (-m, f3x, (3% ; ipi, <p 2 \ x *)> 



where 



(38) 



01 


= 1 - 


to — c 


02 


= 1 - 


m - px 


<Pl 


= 1 - 


m — a\ 


<P2 


= 1 - 


m — a 2 


®m 


(ax 


to) (a 2 \m) (— 




(c m) (pi to) 



Let us first separate the even/odd sectors of J AC as follows: For I — even = 2m , m = 0, 1,2, ... we define 



jpeven _ ^Jpg 



3^2 



-to, -m + \, j + D;l + i+j + D,l-i-2m-D 



and for I = odd = 2to + 1 , to — 0, 1, 2, ... we define 

1 



3^2 — 3-^2 



-, j + D ; 1 + i + j + L>, -i -2m- D 



Using equation (38) in equation (f40D we get 



3^2 



— T 3^2 



-i - i - ^ - A i + \l + D ; 1 - j - -/ - D, ~ 



where 



r^i (-¥ + h\h l )U + D \h l ) 

( ' l + i+j + D\\l)(l-i-l-D\\l) 
(j + D\\l) {-% -j-D\-\l)(i + l + D\- \l) 



ai 



Plugging this into the expression of J z we get 

J AC 



x 3 F 2 [ -ll 



-i -3 



T{-i)T{-j)T{i+j + l + D)T{\) 
\l-D,i + h + D-l- 3 -\l-D, I 



Using the duplication formula for the gamma function 



r(i)r(i + i) 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



even 
AC ' 



2'r(i + ±o 

and rearranging the gamma funtions in convenient Pochhammers' symbols, we arrive at the expression for J A 
i.e., j ac — j ac 

In a completely analogous way, starting from 3 F 2 dd we arrive at J A q . 

Therefore, when equation (|38l ) is applied to our case in J AC with % = 1 it lead us to the following conclusion: The 
I — even sector of J AC yields exactly the y = even sector, J A ^ ven , whereas the I = odd sector of J AC yields exactly 

the y = odd sector, J A Q dd . In order to arrive at these identities, in the intermediate steps of the calculation one needs 
to use the duplication formula for the gamma function. 

Another identity between 3 F 2 hypergeometric functions of unity argument is given by 
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where F p (0; 4, 5) is defined in ( |36[ ) and 

F p (0; 2, 3) 



F p (Q; 4, 5) = F p (0; 2, 3) 



3 f 2 (e - a, f - a, s ; s + 6, s + c 1 1) 

r(a)r( s + 6)ro + c ) 



(46) 



(47) 



Plugging in the parameters of the 3F2 hypergeometric function (see table in section III) into (|36j), the identity ( |46| ) 
above yields 



1, 



1 



Fn = M 3F2 \ 1 + i, 1 — a 1 — D, -— D:- — -I — D,2 — -I — D 



3^2 

where M is a factor given by ratios of gamma functions 

M 



r(!-D)r(i + * + j + i>)r(i 



l-D) 



r(j + D)T(l-U-D)T(2-y-D) 



(48) 



(49) 



If we now redefine the 3F2 hypergeometric function on the right hand side of pq ) to be our new F p cw (0- 7 4, 5) and 
using the fact that for terminating series the following identity is valid 



where 



then we have 



T(s) T(e - c) T(/ - c) F p ncw (0; 4, 5) = T(l - a) T(l - f + b) T(l - e + b) F p (l; 0, 2) 
3 F 2 (l-o ) 1-/ + 6, l-e + 6;2-«-o, 1 — a + 6 1 1) 



F p (l; 0, 2) = 

3^2 = A/" 3-F2 



r(c)r(2-s-a)r(l-a+6) 



(50) 
(51) 



(52) 



We do not need to be overly concerned about the AT factor since both M. and M are ratios of gamma functions 
that at the end can be rearranged conveniently to yield the desired factor present in the J\ c solution. Therefore, 
after some algebraic manipulation of this sort we have J Z AC = J%q. 

In a similar manner, if we interchange parameters a and b in F p ew (0; 4, 5) and proceed as above, we get J Z AC — J\c- 
Lastly, for terminating 3^2 hypergeometric series with parameter c = — m, the following identity is verified 



where 



and 



K(3; 4, 5) 



F p (0; 4, 5) = fif„(3; 4, 5), 

= f iyn r(i- a )r(i-b) 
( j r( e -c)r(/-c)' 

3-^2(1 — a i 1 — b, s ; 1 — a — b + e, 1 — a — b 



/|1) 



T(c) r(l - a - c + e) T(l - a - b + f) 
Substituting the parameters of the hypergeometric function 3^2 in ( |53| ) we get 

3 F* =V 3 F 2 (l-j-D, 1 + ^, \-D; 2 + i + ^l, 2-i-j- \l-2D 



(53) 
(54) 

(55) 
(56) 



where V is a ratio of gamma functions with which we are not going to be concerned about. 

Redefining the RHS hypergeometric function in (|5^) as our new Fp CW (0; 4, 5) and using ( |50| ) we conclude that 
3 F£ = Q 3^, so that, at the end, J\ c = J% c . 

This concludes our proofs of degeneracy in the solution for the Feynman integral. 
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